Statistics of Multiple Sign Changes in a Discrete Non-Markovian Sequence 
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We study analytically the statistics of multiple sign changes in a discrete non-Markovian se- 
quence = + {i = 1,2, ...,n) where cj>i^s are independent and identically distributed 
random variables each drawn from a symmetric and continuous distribution p(0). We show that 
the probability Pm(n) of m sign changes upto n steps is universal, i.e., independent of the dis- 
tribution p{4'). The mean and variance of the number of sign changes are computed exactly for 
all n > 0. We show that the generating function P(p,n) — y^^_n Pm(n)p^ ~ exp[— 0d(p)n] 
for large n where the 'discrete' partial survival exponent 9d(p) is given by a nontrivial formula, 
Od{p) ~ log [sin ""^(-^1 — p^)/\/l — p^] for < p < 1. We also show that in the natural scaling limit 
m — » CX3, n — » oo but keeping x = m/n fixed, Pm{n) '-^ exp[— n<l>(x)] where the large deviation 
function ^{x) is computed. The implications of these results for Ising spin glasses are discussed. 
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The probability Po{T) that a stochastic process il'iT) 
does not cross zero upto time T is a quantity of long 
standing interest to both physicists and mathematicians 
[0,^ and has resurfaced recently with a new name 'per- 
sistence' in the context of nonequilibrium systems [^. 
A lot of recent efforts have been devoted to compute 
Po(T) for stationary Gaussian processes. Such a Gaus- 
sian stationary process (GSP) is completely specified by 
its two point correlation function C{T) = {'i(j{0)'ip{T)) . 
For a wide class of correlation functions, it is known 
that Po{T) exp{—dT) for large T where the persis- 
tence exponent 9 is usually nontrivial, depends on the 
full function C{T) and is calculable exactly only in very 
few cases |^ . A natural generahzation of Pq (T) is Pm (T) , 
the probability of m zero crossings upto time T. The 
mean and the variance of the number of zero crossings 
upto time T of a GSP have been studied before [Q. 
For a smooth GSP where C{T) = 1 - aT^ + . . . for 
small T with a > 0, the mean is given by Rice's for- 
mula {m)/T ~ —C" {Q) /tt and the variance by a 
more complicated formula due to Bendat Recently 
it was shown Q that for a smooth GSP the generating 
function P(p,T) = E"=o^-(V ^ exp[-0(p)T] for 
large T where the 'partial survival' exponent 9{p) varies 
smoothly from 9{Q) = 9 to 9{\) — as p varies contin- 
uously from to 1. In Ref. g, the exponent 9{p) was 
computed exactly for the one dimensional Ginzburg Lan- 
dau model of deterministic coarsening and also approx- 
imately within the independent interval approximation 
for other smooth processes such as the diffusion equa- 
tion. The only process for which a closed form expression 
of 9{p) exists, so far, is the random acceleration prob- 
lem |7|J^ which corresponds to a GSP ||] with correlator 
C(T) = [3exp(-r/2) - exp(-3r/2)]/2. For this prob- 
lem, the exponent 9(jj) for < p < 1 is given by the 
formula, %) = i [l - f sin"! (p/2)] . 

In this Rapid Communication, we study the statistics 
of multiple crossings or sign changes in a discrete se- 



quence -01, "02, ■ ■ "^n as opposed to a continuous process 
^]{T) discussed in the previous paragraph. This study 
is motivated by the recent works on the persistence of 
a discrete sequence [p|-p|. The principal motivation for 
studying the persistence of a discrete sequence is twofold. 
First, in various experiments and numerical simulations 
to measure the persistence Pq (P) of a continuous stochas- 
tic process ip{T), one usually samples the continuous pro- 
cess only at discrete time points separated by a fixed 
window size AT and checks whether the process has re- 
tained the same sign at all these discrete times. Some 
information gets lost due to this discretization since the 
continuous process ip{T) may have crossed and recrossed 
zero in between two successive discrete points. Thus the 
'discrete-time' persistence Po{n) [i.e., the probability that 
the sequence ^'(0), V'(AT), V'(2AP), . . ., ^'(nAP = P) 
have the same sign] is usually greater than the continu- 
ous time persistence Pq (P) . In Ref. |^ , it was shown that 
Po{n) ~ exp[— ^d'T-] where the exponent 9d depends con- 
tinuously on the window size AP. The second motivation 
for studying the persistence of a sequence follows from 
the observation that many processes in nature such 
as weather records are stationary under translations in 
time only by an integer multiple of a basic period (which 
can be chosen to be unity without loss of generality). It 
was shown in Ref. ||l^ that for a wide class of such pro- 
cesses, the continuous time persistence Po(P) is the same 
as the persistence Pq (n) of the corresponding discrete se- 
quence obtained from the measurement of the process 
only at integer times. A natural generalization of Po{n) 
is clearly Pm{n), the probability that there are m sign 
changes along a sequence of size n. 

The exact calculation of Pm{n) for an arbitrary sta- 
tionary sequence seems difficult. It is therefore important 
to find exactly solvable cases. In this paper we present 
exact results for Pm (n) for a specific sequence which was 
introduced in Ref. pO| , 
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where (j>{i)^s are independent and identically distributed 
(i.i.d) random variables, not necessarily Gaussian, each 
drawn from the same symmetric continuous distribution 
p{(t>). The variables ipi^s have only nearest neighbour cor- 
relations. The sequence in Eq. ([^) is stationary but non- 
Markovian since depends not just only on ipi^i but 
on all the preceding members of the sequence . This 
sequence appears as a limiting case of the diffusion equa- 
tion on a hierarchical lattice jlOj . It also appears in the 
one dimensional Ising spin glass problem where tpi rep- 
resents the energy cost to flip the i-th spin In Ref. 
pO| , the persistence Po{n) for this sequence was com- 
puted exactly for all n and remarkably Po{n) was found 
to be universal, i.e., independent of the distribution 
In particular, it was that Poin) ~ exp[— ^^n] for large n 
with 0d — log[7r/2]. The persistence Po{n) was shown to 
be identical to the average fraction of metastable config- 
urations (originally computed in Ref. ||l^) in the corre- 
sponding Ising spin glass chain [0 . 

The purpose of this paper is to show that Pm{n) for 
any m > can also be calculated exactly for the sequence 
in Eq. (|l|) and turns out to be universal. Let us sum- 
marize our main results which are all independent of the 
distrbition p{4>): 

• We show that the mean number of sign changes upto 
n steps (i.e., when the sequence size is n -f 1) is given by 
the exact formula, {m) = n/3 for all n > 0. 

• The variance is given by the formula, cr^ — {m?) — 
{mf = [16n + 3 + 5n,i]/90 for all n > 0. 

• We show that analogus to its continuous counterpart, 
the generating function P{p,n) — X]m=oP™^™("') ^ 
eyL\ji[—9 dij)) n] where the 'discrete partial survival' expo- 
nent 9d{p) is given by the closed form expression 



dd{p) = log 



< p < 1. 



(2) 



This result can be analytically continued to p > 1. 

• We also show that in the limit m — > oo, n — > oo but 
keeping x = m/n fixed, Pmin) ~ exp[— n$(a;)] where 
$(x) is a universal large deviation function that we com- 
pute. 

We start by defining P^rX't'o) to be the joint proba- 
bility that the first member of the sequence in Eq. (|l|) 
±?/)i > and that the sequence undergoes m sign changes 
upto n steps, given the value of (/jq. It is then easy to see 
that they satisfy the following recursion relations 



n;,„+l(0o) = / d0lp(</.l)F+ „((/.!) 

P^,n+i{M = / d(/)ip(0i)F„ (0i), (3) 



where F^ J(j)i) = P^^Mi) + ^m-i,„('/'i) and the ini- 
tial conditions are P;^o((^o) = for m > 0, P^^icjio) = 



P(0i)d0i and Po,o('^o) = 1 - Pofi{M- The generat- 
ing functions P^{p,(t>o) = YlV PmjMo)p"^ then satisfy 
the recursion relations 

/oo 
rf01p((/.l)i^+(p,01) 
-00 

P-+MM= I dcl,,pi^,)F-{p,cj,^), (4) 



where F±(p, 0i) = P^f{p,(t)i) + pP^{p,(t)i) with the 
initial conditions, P^{p,(j)Q) = p((j)i)d(j)i and 

P~{p,(f>o) = 1 — P^{p,<j>o). Further simplification can 
be made by differenting Eq. (||) with respect to (po fol- 
lowed by a change of variable (f)o u = J^" p{(f))d(j) 
and then using the symmetry p{(j3) ~ p{—(p). Writing 
Pn{Pi4>o) — pti'^) (suppressing the p dependence for 
convenience) we find two coupled non-local recursion re- 
lations 



pp 



(5) 



with the initial conditions p^{u) = 1/2 ± u and the 
boundary conditions p*(=Fl/2) = 0. Note that the ex- 
plicit dependence on the distribution p{(j)) disappears in 
Eq. (||). As a result, all further quantities computed from 
these recursion relations will be independent of p{4i) pro- 
vided p('/>) is symmetric and continuous. 

In principle, one can solve the recursion relations in 
Eq. (^ by the generating function method. However to 
calculate the mean and the variance of zero crossings, it 
is simpler to directly analyze Eq. (|^). Let £^^(0o) — 
12m=o Pm n('^o) bc the cxpccted number of sign changes 
upto n steps with the first member ipi positive (or neg- 
ative) and given (pQ. Let us write E^{(j)o) — e^{u) after 
making the change of variable (j)Q u. The average num- 
ber of crossings is then given by (m) = J^^[E^{(f>a) + 

1 /2 

/2[^n (") + {u)]du. Differenti- 
ating Eqs. 
we get 



once with respect to p and putting p = 1, 



de 



n+l 



du 



= ± [e+(-M) + e„(- 



u ± — , 
2' 



(6) 



with the initial conditions e^(u) = and the bound- 
ary conditions e^(=Fl/2) = 0. These recursion rela- 
tions can be solved exactly and one gets e^^^ — (1/2 ± 
M)(n-4)/12±u3/3-Hu2/2±u/2-h 1/6. Hence e„(u) = 
e+(u) + e^(u) = + ?i/3 - 1/12. This then gives the 

1/2 

exact result (m) = J_-^^^en{u)du — n/3 for all n > 0, 
independent of p{(t)). 

The variance cr^ = (m^) — (m,)^ can be computed in 
a similar way by differentiating Eq. (^ twice with re- 
spect to p, putting p = 1 and solving the resulting recur- 
sion relations. The functions Gfi{(f>o) = X]m=o ™('^ ~ 



2 



i)Pm.niM = 5,t(^) satisfy the following inhomogeneous 
recursion relations for n > 0, 



dg„ 



du 



(7) 



with the initial conditions g^{u) = 
boundary conditions g^{^\/2) = 0. 



gfiu) = and 
Using the known 



values of e„(w) one can again solve Eq. M) exphcitly. 

1 /2 

This finally gives (m(m-l)} = J_y^[g+ {u)+g:;r^ {u)]du = 
(10^2 - 14n + 3)/90 for aU n>2. Using (m) = n/3, one 
gets cr^ = [I6n + 3 + (5„.i]/90 for all n > 0, again inde- 
pendent of 

We now turn to the calculation of the partial survival 
exponent. We expect that for large n, pt{u) ~ f^{u) 
where A = exp[0d(p)]. Substituting this asymptotic form 
in Eq. we get the non-local eigenvalue equation 



du 



±A [f+{-u) + r{-u)\ 



(8) 



subject to the two boundary conditions, /+(— 1/2) = 
and /~(l/2) — 0. Diagonalizing Eq. (Q) and solving the 
resulting non-local equations we get the most general so- 
lutions of Eq. d) 



f^iu) = cos(^u) + sin(^M), 



(9) 



where fj, = X\/l — and the four constants and 
can be written in terms of only two unknown constants 
a and h via the relations, A'^ = ap, A^ = b^l — — a, 
_B+ — bp and = a^l — p'^—b. The solution in Eq. (||) 
must satisfy the two boundary conditions /^(=Fl/2) = 
which gives two homogeneous linear equations for the 
unknown constants a and b. Eliminating a and b from 
these two equations one gets /i = sin~^[yA^jP] and 
hence A = sm^^[^yl —p^\l^\ —p^. Using the relation 
Bdip) — log A, we obtain the result in Eq. (||), once again 
independent of Note that for p = 0, reduces 

to the usual discrete persistence exponent 9^ = log(7r/2). 

The expression for 9 dip) in Eq. is valid in the 
range < p < 1. However, in principle, one can define 
the generating function P{p,n) = Yl'^^o Pra{n)p'^ even 
for p > 1. Then for p > 1 one expects -P(p, n) to diverge 
as n — > 00 indicating 9d{p) becomes negative for p > 1. 
Indeed one can easily get the result for p > 1 by analyt- 
ically continuing the expression in Eq. (|^) to the range 
p > 1 and this gives 



%{p) = log 



log (p + - 1 



(10) 



Thus 9d{p) tends to 
log[log(2p)/p] as p ^ oo. 



-oo rather slowly as 9d{p) 



We next analyse the distribution Pm{n) in an inter- 
esting scaling limit. Since the average number of cross- 
ings scale linearly with the size n as (m) — n/i, a nat- 
ural scaling limit is when m — > oo, n — * cxd but keep- 
ing the ratio x — m/n fixed. In this limit, we show 
that Pm{n) ~ exp[—n<I>(a;)] where $(a;) is a large devia- 
tion function which is universal, i.e., independent of the 
distribution /3(0). Large deviation functions associated 
with different physical observables have appeared before 
in the context of various nonequilibrium systems [ p^ . 
The present model provides an example where there is 
a large deviation function associated with the number of 
zero crossings that can be computed explicitly. Indeed, 
substituting the ansatz Pm{n) ^ exp[— n$(a;)] in the gen- 
erating function one gets P{p,n) ^ Ylm=o Prn{n)p"^ ^ 
da; exp[— n{$(a;) — a:logp}]. In the large n limit 
the integral can be evaluated by the steepest descent 
method and one gets P(p, n) ~ cxp[— nG'(p)] where 
G(p) — mina;{$(x) — xlogp}. On the other hand, by def- 
inition, P{p,n) ~ eiq)[—9d{p)'n] for large n. This estab- 
lishes the relation, mmx{^{x) — xlogp} — 9d{p)- Thus 
Od{p) is just the Legendre transform of ^(x). Inverting 
this Legendre transform we get 



^{x) = maxp [xlogp -I- 9d{p)] 



(11) 



where 9d{p) is given exactly by Eqs. (||) and (p^ ) in 
the range p > 0. Note that determining ^{x) from Eq. 
( |ll| ) requires a knowledge of 9d{p) not just in the range 
< p < 1 but also for p > 1. Thus we will need both the 
formulas in Eqs. (0) and (|o|). 
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FIG. 1. The large deviation function plotted against 

X. The solid line represents the function obtained using Math- 
ematica. The dotted line represents the analytical asymptotic 
form ^{x) = — log(l — a;) — — a;) log 2 in the limit x ^ 1. 

The function ^{x) log(7r/2) = 0.451583 ... as a; ^ 0. 
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We have obtained $(a;) from Eq. ( pi] ) using Mathe- 
matica and is displayed in Fig. 1. Since the number of 
crossings m < n, the allowed range of a; is < x < 1. 
One can analytically determine the behavior of $(a;) in 
the three limits x ^ 0, x ^ 1 and x 1/3. First con- 
sider the limit x — > 0. This corresponds to p limit 
of Od{p)- Expanding Eq. (||) for small p, we get Od{p) ~ 
log(7r/2 —p). Substituting this in Eq. ( pi] ) and maximiz- 
ing with respect to p gives $(a;) ~ log(7r/2) -f a;log(a;) as 
a; — > 0. Next consider the opposite limit when a; — > 1. 
This limit correspond to 0dip) in the limit p —> oo. Hence 
we need to now use the analytically continued formula 
in Eq. (|o|). Expanding Eq. (^ for large p, we get 
Od{p) « log[log(2p)/p] to leading order. Substituting this 
asymptotic form in Eq. (|l^) and maximizing with re- 
spect to p we get $(x) « — log(l — a;) — 1 + (1 — x) log 2 
as X ^ 1. In Fig. 1, this asymptotic form is shown by 
the dotted line to which $(x) approaches rather quickly 
as X ^ 1. 

The most interesting limit, however, is when x —^ 1/3, 
i.e. m — > {m). This limit in x corresponds to p ^ 1 limit 
of Od{p). It is easy to see that both the limits p ^ 1~ 
and p — > 1"*" yield the same result. Let us consider the 
case when p ~ 1 ~ e where e ^ 0. Expanding Eq. (H) in 
powers of e, we get 0d{p) « e/3 -I- 7e^/90 + O(e^). Substi- 
tuting this in Eq. (|ll|) and maximizing with respect to 
p = 1 — e, we get as x — > 1/3, 

(12) 



$(x) « - 



This limiting form can also be derived independently 
from a central limit theorem. To see this we write the 
number of sign changes m as the sum m — X)r=i 
with Wi = 1 — 0{'ipi'4>i+i) and 0{x) is the Heaviside step 
fimction. Thus m — (m) = J27=i('^i ~ i'^i))- Clearly 
in the limit m — > (m), the variables {wi — (wi)) be- 
come only weakly correlated. Then in the limit when 
n is much larger than the correlation length between 
these variables one expects the central limit theorem to 
hold predicting a Gaussian distribution for m, P„i{n) ~ 
exp[— (m — (m))^/2cr^]. Using the already derived re- 
sults (to) = n/3 and cr^ « 8n/45 for large n, we find 
Pm{n) ~ exp[— 45ri(x — 1/3)^/16] thus yielding the same 
$(x) as in Eq. (p^). Thus this limit provides an indipen- 
dent check of our results for the mean and the variance. 
The three limiting behaviors of $(x) are summmarized 
as follows 

log(7r/2) -I- xlogx X ^ 0, 

"3) X — > 1/3, 

-log(l - x) - 1 + (1 - x)log2 X ^ 1. 

(13) 

We conclude with a discussion of the implications of 
our results for an Ising spin glass chain described by the 



$(x) 



Hamiltonian, H = Ji^i^iSiSi+i with Si = ±1 and 
the bonds J^^i+i's are i.i.d random variables each drawn 
from the same symmetric and continuous distribution. A 
spin will be called metastable if the cost of energy to flip 
it under zero temperature Glauber dynamics is positive, 
i.e. /\Ei = 2s.i[Ji-i^iSi-i + Ji^i+iSi+i] > 0. A given spin 
configuration (with fixed J's) consists of alternate do- 
mains of metastable and non- metastable spins. A natural 
question is what is the average (over disorder) probabil- 
ity P(to, n) that there are to such domains in a chain of 
length n. Defining the new variables = 2Ji^i+iSiSi+i 
which are also i.i.d variables, we see that the energy costs 
l^Ei = form exactly the sequence studied in this 

paper. The average domain number probability P(m, n) 
is then identical to the probability of having 2to sign 
changes in the sequence {AEi} upto n steps, i.e., P2m{n) 
that has been computed exactly in this paper. Clearly 
for TO = 0, P{0,n) = Poin) is just the fraction of fully 
metastable configurations (out of the 2" configurations) 
at zero temperature and is the same as the persistence of 
the sequence {AEi} (lO). Note that one can easily gen- 
eralize this question of the domain number probability 
to higher dimensions as well. The study of the statistics 
of the domains of metsatable spins in higher dimensional 
spin glasses may provide interesting insights into the na- 
ture of the low temperature phase. 

I thank D. Dhar and A.J. Bray for useful discussions. 
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